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Abstract 

The PARC2D code has been selected to analyze the 
flowfields of a representative hypersonic scramjet nozzle 
over a range of flight conditions from Mach 3 to 20. The 
flowfields, wall pressures, wall skin friction values, heat 
transfer values and overall nozzle performance are pre- 
sented. 

Nomenclature 
a = sound speed 

Cf = skin friction coefficient = C/ = 

Ch = heat transfer coefficient = Ck = 

C p = specific heat 

C v = velocity coefficient 

r] = nozzle efficiency = r? = 

7 = ratio of specific heats 

h = enthalpy 
M = Mach number 
[i = viscosity 
n = unit normal 
P w = wall pressure, p w /p c 
P T — Prandlt number 
R e = Reynolds number 
p = density 
T = temperature 
X = streamwise distance 

Subscripts 
c = combustor exit 
e = nozzle exit 
1 = laminar 

o = freestream condition 
s — constant entropy 
t = total condition, turbulent 
w = wall 


* Supervisor, Hypersonics Analysis Section, Senior Mem- 
ber AIAA 

** Senior Research Engineer, Hypersonics Analysis Sec- 
tion, Member AIAA 

*** Research Engineer, Hypersonics Analysis Section, 
Member AIAA 


Introduction 

Accurate evaluation of nozzle performance is essen- 
tial in hypersonic propulsion due to the high sensitivity 
of net thrust to nozzle performance. The nozzle may be 
highly integrated with the airframe so that the upper wall 
of the nozzle is actually the aft end of the airframe. The 
lower wall can be cut shorter than the upper wall to save 
weight. Examples of these types of nozzles have been given 
by Lewis et. al. 1 , and Lilly 2,3 . 

An existing computer code, PARC2D 4 , is selected 
for this study. The code solves the full two-dimensional 
Reynolds averaged Navier-Stokes equations in strong con- 
servation form with the Beam and Warming approximate 
factorization algorithm. The implicit scheme uses cen- 
tral differencing for a curvilinear set of coordinates. The 
code was originally developed as AIR3D by Pulliam and 
Steger 5 , and Pulliam 6 later added the Jameson 7 artificial 
dissipation and called the code ARC3D. Cooper 4 adapted 
the code for internal propulsion and called the code PARC. 
The current study was conducted in 2-D and could have 
been conducted in 3-D with an increase in computer re- 
sources required. Previous configurations examined at 
NASA Lewis by Reddy and Harloff 8 include: 1) a 2-D 
and a 3-D hypersonic nozzle at Mo = 5, comparisons were 
made with computations of Spradley 9 and others; 2) a 2-D 
ramp at Mo = 14.1 with ramp angle of 18°, Holden 10 ; and 
3) a 3-D corner flow configuration at Mo = 2.94 of Oskam 
et. al. 11 . The study by Reddy and Harloff 8 provided con- 
fidence in the PARC codes for a variety of complex shapes 
and high speeds. Also Cooper 12 , obtained reasonable flow- 
fields with PARC2D for trumpet and conical nozzles. 

The goal of the current study is to examine the flow- 
fields and overall nozzle performance over a wide range 
of flight conditions. There are no available experimental 
data, with external flow, to compare with the present com- 
putational results. 

Nozzle Geometry 

Figure 1 shows the nozzle geometry used in this study. 
The initial divergence angle is about 37°. The geometry 
was obtained from McFarland 13 . 
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Freestream and Combustor Exit Properties 

The ratio of specific heats, gamma, and the gas 
constant of the combustor exit are considered constant 
throughout the entire flowfield and the values are shown 
in Table 1. The static nozzle entrance Reynolds numbers, 
based on nozzle entrance height, are such that turbulent 
flow is expected for the entire Mach number range studied 
and turbulent flow is assumed. The freestream and com- 
bustor exit properties were obtained from McFarland 13 . 

Grid 


The grid used for the computations was 199 x 129, 
uniform in the X direction, see Figure 2. Grid clustering 
near the upper and lower walls is needed to resolve the 
boundary layers. Also grid clustering downstream of the 
cowl is needed to resolve the shear layer. The present study 
did not adapt the grid to the shear layer location, although 
this has been accomplished by the authors in previous un- 
published studies. The solutions obtained with adapted 
grids were not significantly different. 

The boundary conditions are uniform inflow and ex- 
trapolated outflow properties for supersonic flow. The su- 
personic/hypersonic inflow boundary above the exit was 
assumed to be uniform. For subsonic outflow the pressure 
is imposed by the sublayer approximation. These bound- 
ary conditions are ideal and more realistic inflow profiles 
can be accounted for. The wall temperature is assumed to 
be constant at 1750 deg R. 

Mach 3 Flowfield 

The velocity vectors at Mo = 3 are shown in Figure 
3. The flow is overexpanded and a shock propagates from 
the cowl lip to the upper nozzle wall. The exhaust flows 
along the upper nozzle wall similar to a wall jet. A slight 
separation occurs on the upper wall ahead of the shock. 
This is difficult to see in this figure, and it will be further 
discussed below. The Mach number contours are shown in 
Figure 4. Due to the overexpanded nature of the flow at 
the end of the lower cowl, the shear layer turns upward. 
The shear-layer deflection and pressure imbalance cause a 
shock to form above the shear layer. 

The total pressure contours are shown in Figure 5. 
The gradient in nozzle and freestream total pressure across 
the shear layer is illustrated. As mentioned previously, the 
shear layer deflects upward due to the overexpansion of the 
nozzle flow at the cowl exit. 

The divergence of velocity provides a scalar represen- 
tation of density gradients such that compressive and ex- 
pansive regions of the flow can be visualized. Compres- 
sion regions are dark and expansion regions are light. The 
shock from the end of the cowl to the upper nozzle wall is 
illustrated in Figure 6 as a dark region. The shock deflec- 
tion angle is about 39° and the Mach number in front of 
the oblique shock is about 3.5. The shock reflects off the 
boundary layer and then propagates in the downstream 
direction, i.e. left to right. The boundary layer separates 


on the upper nozzle wall ahead of the incident shock. As 
the shock interacts with the wall boundary layer the shock 
initially reflects at a sharp angle. As the flow moves over 
the separation bubble, the flow turnes less, as shown by 
the decrease in magnitude of the dark-colored compres- 
sion region. At the aft end of the separation bubble, the 
flow is again turned at a sharp angle into the reflected 
shock downstream of the separation. In the lower region 
of the flowfield, below the cowl, a shock emanating from 
the leading edge of the cowl is observed. This shock is 
caused by flow deflection due to the boundary layer devel- 
opment along the lower cowl wall. The expanding region of 
the flow downstream of the cowl is observed in the Figure 
6 as a light region. 

The wall pressures for the three nozzle wall surfaces, 
see Figure 7-a, initially decrease with increasing X for both 
upper and lower nozzle walls. The interior nozzle shock 
hits the upper nozzle wall at about 2.8 feet as indicated by 
in the marked pressure rise at this location. This feature 
has been shown in previous figures. The pressure below 
the cowl is larger than the pressure in the nozzle and this 
pressure difference at the end of the cowl causes the shock 
to form in the nozzle. The wall skin friction values are 
shown in Figure 7-b for the three surfaces. The skin fric- 
tion on the upper nozzle wall goes negative at a location of 
about 2.8 ft indicating a small region of reverse flow. The 
skin friction of the upper cowl wall is slightly higher than 
that of the upper wall, and is also higher than the skin 
friction of the lower cowl wall. The heat transfer for the 
upper and lower nozzle walls are similar in magnitude over 
the length of the cowl, see Figure 7-c. A slight decrease, 
then a rise in the heat transfer coefficient is observed near 
X = 2.8 ft. where the boundary layer flow separates. In 
the nozzle, the heat transfer is positive indicating hot gas 
and cold wall. Below the cowl the heat transfer coefficient 
is negative indicating cold gas and hot wall. 

Mach 6 Flowfield 

The velocity vectors at Mach 6 are shown in Figure 8. 
The flow is slightly overexpanded at the cowl exit and the 
flow expands behind the cowl by means of an expansion 
fan. The wake is deflected upward as is indicated by the 
velocity vector angles and the position of the shear layer. 
Reverse flow at the upper boundary is not indicated. The 
Mach number contours for the flowfield, shown in Figure 
9, also illustrates the location of the shear layer. The total 
pressure contours, see Figure 10, show the upward deflec- 
tion of the shear layer. The divergence of velocity are 
shown in Figure 11. The dark region near the upper wall 
differs from that of the Mach 3 flowfield. This compression 
region is believed to be associated with a coalescence of 
the shock from the lower cowl wall and a standing oblique 
shock in the nozzle exit. The wall pressures are shown in 
Figure 12-a and the overexpansion is illustrated. The wall 
skin friction and heat transfer values are shown in Figure 
12-b and 12-c respectively. 
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Mach 10 Flowfield 


Summary 


The exhaust flow is underexpanded at Mach 10 as in- 
dicated by the downward deflection of the nozzle flow, see 
Figure 13. There are no indications of shocks inside the 
nozzle. A shock propagates from the lower cowl lip down- 
ward, due to the underexpansion of the nozzle flow, see 
the Mach contours in Figure 14. Figure 15 shows the total 
pressure contours at Mach 10. The downward deflection of 
the shear layer, is illustrated. The shock from the end of 
the cowl to the region below the nozzle is illustrated as a 
dark region below in Figure 16. The shock deflection angle 
is about 18°. The pressure is constant on the lower wall, 
see Figure 17-a. The C / and C* wall values are shown in 
Figure 17-b and 17-c respectively. 

Mach 20 Flowfield 

The velocity vectors at Mach 20 are deflected down- 
ward aft of the cowl due to the underexpansion, see Figure 

18. A shock propagates downward from the lower cowl. 
The lower cowl leading edge shock merges with the trail- 
ing edge shock about 26% of a cowl length downstream of 
the cowl. 

The flow along the curved upper nozzle surface turns 
and compresses, like flow along a curved ramp, and a 
curved shock forms inside the nozzle. An upwardly de- 
flected shock behind the upper nozzle wall is evident in 
both Figure 18 and the Mach number contours of Figure 

19. The total pressure contours in Figure 20 show the 
downward deflection of the shear layer. The shock from 
the lower side of the cowl to the region below the nozzle is 
illustrated in the divergence of velocity contours, see Fig- 
ure 21. The shock from the leading edge of the cowl merges 
with the shock from the aft end of the cowl about 18% of 
the cowl length downstream of cowl. 

An interior shock along the upper wall is evident in 
Figure 21. This interior shock is typical at overspeeded 
conditions 14 , i.e. when the entrance Mach number exceeds 
the design value. Notice that the upper wall pressure does 
not follow the lower nozzle wall pressure but instead in- 
creases the full length of the nozzle. The nozzle flow un- 
derexpansion is evident in Figure 22-a. The C f and C p 
values are shown in Figure 22-b and 22-c respectively. 


Nozzle Performance 

The adiabatic nozzle efficiency, * 7 , is the ratio of the 
nonisentropic and isentropic enthalpies between the com- 
bustor exit, c, and the nozzle exit, e. The velocity coeffi- 
cient, C v , is the ratio of the actual to the 1-D isentropic 
velocity and is the square root of the nozzle efficiency. The 
velocity coefficient increases with Mach number up to Mo 
of 10 and then decreases. For example, C v is 0.93 at Mo 
of 3, 0.98 at Mo of 6, 1.0 at Mo of 10 and 0.94 at Mq of 

20. These values might differ with chemistry accounted 
for. This will be checked in future research. 


Flowfields and wall values of P^, Cj , and C h have 
been computed for a scramjet nozzle from Mach 3 to 
20. Overall nozzle performance has also been computed. 
Matching the external pressure at the cowl exit is necessary 
to avoid waves in the nozzle. At Mach 3 and 6 the nozzle 
flow is overexpanded and shocks are inside the nozzle and 
the shear layer is deflected upward. At Mach 10 and 20 
the nozzle flow is underexpanded and shocks are outside 
the nozzle and the shear layer is deflected downward. At 
Mach 20 a shock forms inside the nozzle along the upper 
wall. 
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Table 1. Freestream and Combustor Exit Properties. 


M 0 

T 0 

R 

p 0 

PSIA 

M c 

T c 

R 

p c 

PSIA 

F- 

R 

-LB/R-I 

3 

390 

1.7 

1.0 

4075 

11.9 

1.25 

62.7 

6 

393 

.4 

1.5 

4220 

10.4 

1.19 

63.4 

10 

409 

.15 

1.9 

5403 

75.0 

1.16 

66.0 

20 

468 

.03 

4.6 

5119 

25.9 

1.18 

88.0 


NOZZLE DEFINITION 

Y = A1 + A2 X + A3 X**2 + A4 X **3 

A1= .02304 
A2= .75459 
A3= -.12370 



Figure 1. Nozzle Definition. 
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gure 3. Mach 3, Velocity Vectors. 


A 0.00Q0 
B 0.9320 
C 0.8691 
□ 1.2961 
E 1.7282 
F 2.1602 
G 2.5922 
H 3 .024*3 
1 3 • *563 
J 3.8883 



A 0.0073 
B 0.3577 
C 0.7082 

0 1.0587 
E l.*092 
F 1.7596 
G 2.1101 
H 2. *606 

1 2.8110 
J 3.1615 



Figure 5. Mach 3, Total Pressure Contours. 
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Figure 6. Mach 3, Divergence of Velocity. 
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AXIAL DISTANCE 
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Figure 7. Mach 3, P w , Cj , C h vs. X. 
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Figure 8, Mach 6, Velocity Vectors. 


A 0.0000 
B 0.7112 
C 1.4223 

0 2.1335 
E 2.8446 
E 3.5558 
G 4.2670 
H 4.9781 

1 5.6893 
J 6.4004 



Figure 9. Mach 6, Mach Number Contours. 
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A 0*0065 
B 15*0583 
C 30. mo 
D **5.163? 
E 60.216** 
F 75.2691 
G 90.3219 
H 105.37H6 
1 1 20. **273 
J 1 35.9801 



Figure 10. Mach 6, Total Pressure Contours. 
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Figure 12. Mach 6, P w , Cf, Ch vs. X. 



Figure 13. Mach 10, Velocity Vectors. 


A 0.0000 
B 1.0177 
C 2.0353 

0 3.0530 
E *.0706 
F 5.0083 
G 6.1060 
H 7.1236 

1 8.1*13 
J 9.1589 



Figure 14. Mach 10, Mach Number Contours. 




A 0.0003 
8 2*0. *279 
C *80.8557 

0 ?21 ,2832 
E 961.7112 
F 1202.1389 
G 1**2.566? 
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1 1923. *221 
J 2163.8*99 
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Figure 17. Mach 10, P w , Cf, Ch vs. X. 


Figure 18. Mach 20, Velocity Vectors. 


A 0.0000 
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Figure 19. Mach 20, Mach Number Contours. 
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Figure 20. Mach 20, Total Pressure Contours 
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Figure 21. Mach 20, Divergence of Velocity. 
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